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ABSTRACT 

Gamma-ray bursts (GRBs) and following afterglows are considered to 
be produced by dissipation of kinetic energy of a relativistic fireball and 
radiation process is widely believed as synchrotron radiation or inverse Compton 
scattering of electrons. We argue that the transfer of kinetic energy of ejecta 
into electrons may be inefficient process and hence the total energy released by a 
GRB event is much larger than that emitted in soft gamma-rays, by a factor of 
~ (m p /m e ). We show that, in this case, very strong emission of TeV gamma-rays 
is possible due to synchrotron radiation of protons accelerated up to ~ 10 21 eV, 
which are trapped in the magnetic field of afterglow shock and radiate their 
energy on an observational time scale of ~ day. This suggests a possibility 
that GRBs are most energetic in TeV range and such TeV gamma-rays may 
be detectable from GRBs even at cosmological distances, i.e., z <~ 1, by 
currently working ground-based telescopes. Furthermore, this model gives a 
quantitative explanation for the famous long-duration GeV photons detected 
from GRB940217. If TeV gamma-ray emission which is much more energetic 
than GRB photons is detected, it provides a strong evidence for acceleration of 
protons up to ~ 10 21 eV. 



Subject headings: acceleration of particles — gamma rays: bursts — gamma rays: 
theory 



1. Introduction 

Gamma-ray bursts (GRBs) are widely believed as dissipation of kinetic energy of 
relativistic motion produced by an expanding fireball with a Lorentz factor of ~ 10 2 -10 3 
(see e.g., Piran 1994 for a review). The recently discovered afterglows following GRBs are 
also considered as similar phenomena, which are dissipation in the external shock generated 
by the collision with interstellar matter (Paczyhski & Rhoads 1993; Katz 1994; Meszaros 
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& Rees 1997; Vietri 1997a). The cosmological origin of GRBs is now almost confirmed by 
the discovery of metal absorption lines at z— 0.835 for the optical afterglow of GRB970508 
(Metzger et al. 1997), and some of observations for X-ray, optical, and radio afterglows are 
in rough agreement with the prediction of the cosmological fireball model (Wijers, Rees, 
Meszaros 1997; Waxman 1997a, b; Vietri 1997b). However, there is large variation in the 
afterglow response of GRBs (e.g., Groot et al. 1998), and it is not yet clear whether the 
simple afterglow model is applicable for all of GRBs. 

There are two important, but highly uncertain parameters in such theoretical models 
of GRBs and afterglows: the degree of equipartition between the internal energy of shock 
heated matter and magnetic fields and between protons and electrons (£ e ). In most 
of publications which calculated model predictions of GRBs or afterglows, these two 
parameters are assumed to be of order unity, and the radiation process is considered as 
electron synchrotron (or inverse Compton scattering). In this case the efficiency of energy 
release in GRBs or afterglows compared to the total energy of a GRB event (E) is of 
order unity. However, currently there is no clear evidence for efficient energy transfer into 
electrons and magnetic fileds, although some of observational data are consistent with 
£ e ~ 1 (Waxman 1997a). If the energy transfer from protons into electrons is inefficient, 
energy stored in electrons is only a fraction of £ e ~ (m e /m p ) of the total fireball energy 
and hence about 2000 times larger energy must be released as kinetic energy of relativistic 
ejecta than the observed energy emitted as GRB photons. The GRB photon energy is 
~ 10 51 erg if the radiation is isotropic and the redshift of most distant GRBs, ^ max ~ I. 
Then the total energy E may be uncomfortably large because most of GRB models are 
based on gravitational collapses of massive stars in which available energy is ~ 3 x 10 53 
erg and most of this energy will be lost as neutrinos. However, the theoretical estimate of 
merger rate of binary neutron stars is about 10 2-3 times higher than the observed GRB rate 
(Lipunov et al. 1995; Totani 1997,1998), which suggests that GRB is strongly beamed if 
GRBs are associated to merger of binary neutron stars (Blinnikov et al. 1984). If GRBs are 
actually beamed with such a strong beaming factor, the above constraint of energy budget 
becomes much weaker. Much more energetic models of GRBs have also been proposed such 
as the microquasar model, in which total energy of ~ 10 54 erg can be supplied to a fireball 
(Paczyhski 1998). 

In this letter we argue that time scale of energy transfer into electrons by the Coulomb 
interactions is much larger than the expansion time of external shock, while magnetic field 
may achieve the equipartition with protons in the shock heated matter. We then show 
that, as a consequence of this scenario, a very strong TeV emission is expected during a 
few days after GRBs by synchrotron radiation of 10 20 eV protons and it may be detectable 
by current ground-based telescopes even from a GRB at cosmological distances, in spite of 
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significant attenuation due to e creation with intergalactic infrared photons. Synchrotron 
emission of protons of ~ 10 20 eV from GRBs was first considered by Vietri (1997c), and 
Bottcher & Dermer (1998) extended the analysis to emission from afterglows. Both papers 
concluded that TeV gamma-rays are detectable only for nearby GRBs (z < 0.1), assuming 
that total fireball energy is of the same order with that of GRB photons. The natural units 
with c = h — 1 is used in this letter. 



2. Efficiency of Energy Transfer into Electrons and Magnetic Fileds 

The evolution of external shock is described by bE = 167rnm p r 3 7 2 /17, where E is the 
total energy released in an opening angle of AQ, b (= 47r / AO) a beaming factor, and 7 the 
Lorentz factor of the shock heated matter (Blandford & McKee 1976). The location of the 
shock, r, is measured in the laboratory frame and n is the (unshocked) interstellar matter 
density. Initially the kinetic energy stored in electrons is only a fraction of (m e /m p ) of the 
total energy, and much greater energy of protons must be efficiently transferred into electrons 
by some interactions in the shock heated matter in order to achieve energy equipartition 
between electrons and protons. However, relative importance of the Coulomb interaction 
becomes smaller with increasing energy of particles, and the particle energy in relativistic 
shocks is much greater than that in non-relativistic shocks such as supernova remnants. The 
time scale of energy transfer in relativistic plasma is difficult to estimate accurately, but a 
rough estimate is given by r ep ~ (nVt) -1 , where n' = A^n is the proton number density 
of the shocked matter measured in the shocked-shell frame and a t = 4irL e (e 2 /m e ^) 2 is the 
transport cross section for electron-proton collisions. The Coulomb logarithm is given by 
L e = ln(am e 7), where a = (m e 7 /Ann'e 2 ) 1 / 2 is the Debye length (e.g., Lifshitz & Pitaevskii 
1981). This time scale should be compared to the expansion time measured in the shell 
frame, r/7, and we find T ep /{r/^) = 1.1 x 10 4 6 _1 / 3 -E 51 1 ^ 3 n ] ^ 2 ^ 3 7 8 / 3 , where E = 10 51 -E 5 i erg 
and n — ri\ cm" 3 . Hence energy transfer through the Coulomb interaction is likely to be 
inefficient. 

On the other hand, magnetic field in the shocked matter may be in equipartition 
with the random motion energy of protons which is directly converted from the kinetic 
energy of a fireball. Recall that, in the well-known equation of magnetohydrodynamics, 
the time evolution of magnetic field, <9B/<9t, is governed by the diffusion term and the 
source term, rot[u x B], where u is the velocity field of fluid. In afterglow shocks there will 
be turbulent motion of u 2 ~ 1 and a coherent length scale will be smaller than the shell 
thickness of the shocked matter measured in the shell frame, ~ r/7. This suggests that 
|rot[u x B]| > {pfjr)B and hence the growth time scale of magnetic field may be < r/7. 
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Since the expansion time is also ~ r/7 in the shell frame, it is possible that equipartition 
between magnetic field and protons is achieved while electrons carry much smaller energy. 
Although the above argument is quite rough and some unknown processes in relativistic 
matter may allow electrons to be in equipartition with protons, it seems rather reasonable 
to consider the case of £b ~ 1 and £ e ~ (m e /m p ) < 1. In order to investigate such an 
energetic model of GRBs, we use E = 10 52 E 52 erg and b = 2006 2 oo as typical values, with 
which E is about 2000 times larger than the energy emitted in GRB photons, i.e., ~ 10 51 6 _1 
erg when 



3. Proton Synchrotron in Afterglow Shock 

If the energy transfer from protons into electrons is inefficient but magnetic field is 
nearly in equipartition, the synchrotron radiation of protons becomes relatively important. 
The energy density of shocked matter is given by A^ 2 nm p in the shell frame and magnetic 
field can be written as B = (32Tc^ B / y 2 nm p ) 1 ^ 2 . It has been considered that protons 
may be accelerated up to ~ 10 20 eV in GRBs because the physical quantities of GRBs 
allow acceleration of protons to such high energies and observed flux of highest energy 
cosmic rays is consistent with the GRB occurrence rate provided that such protons carry 
roughly the same amount of energy with GRBs (Waxman 1995; Vietri 1995). We assume 
that the shock acceleration time is given by ^r^, where = m p , y p /(eB) is the Larmor 
radius, 7 P is the proton Lorentz factor in the shell frame, and r] is a parameter of order 
unity. The maximum energy obtained in the external shock is given by the equation 
V r L = ^lli and we find 4.21 x 10 21 ?7~ 1 ^]/ 2 n^ 6 7ioo&2oo-^52 3 e V in the observer's frame, 
where 7100 = 70/IOO and 70 is the initial fireball Lorentz factor. On the other hand, the 
maximum energy is also constrained by synchrotron cooling. The cooling time at the shell 
frame is t syn = 3mp/(4<T T m 2 f/ mag 7 P ), where t/ mag is the energy density of magnetic field. 
From the equation 7777, = t syn we then find the maximum energy in the observer's frame as 
3.27 x 10 21 ry _1 / 2 ^ B 1//4 n 1 1(/4 7ioo e V, which does not depend on the total fireball energy of 
GRBs. Therefore protons may by accelerated up to io 21-22 eV for 70 = 100-1000, which 
is about one order of magnitude greater than the estimate obtained with bE ~ 10 51 erg 
(Waxman 1995; Vietri 1995). 

The protons accelerated up to ~ 10 21 eV will radiate their energy by synchrotron 
radiation in magnetic fields of afterglow shock. The synchrotron photon energy in 
the observer's frame is given by e 7 = 77 2 e-B/m p = 2.8£ 2 2 i£i/ 2n i^ 2 TeV , where 
£ p ,2i = ' m p77p/(10 21 eV) is proton energy in the observer's frame. Therefore synchrotron 
photon energy is related to proton energy independently of time and will extend to > 
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TeV. Next let us check whether such protons can be confined in the afterglow shock. 
The ratio of Larmor radius of protons to the restframe shell thickness r/7 is found as 
fL/ir/l) = 0.26 G 3/4n r 3/8 ^oo /8 - E 52 3/8 ^TeV t day) where e 1 = e 7)Te v TeV and the observation 
time is tday day = r/(2 7 2 ) or 7 = 15 tda/^oo^^ 8 ^ 1 ^ 8 • The transverse width of beamed 
shell is larger than the shell thickness unless the beaming factor is extremely large as 
b > 7.21 x 10 3 t d " a y£ , 52 3 n^ 1//3 . Therefore the protons which correspond to synchrotron 
photons of < TeV can be trapped within the magnetic field of afterglow shock on a time 
scale of ~ day. The cooling time observed on the Earth is related to the restframe cooling 
time as t syn , bs = t syQ /(2-f), and we find t syn ,obs = 1.3 i^n^e'^t^b^E^ day. 
If the spectrum of accelerated protons is that of the standard shock acceleration theory, 
i.e., dNp/d'jp oc 7"° with a ~ 2, luminosity of synchrotron radiation per decade of photon 
energy, L(e^) = e 1 dL/de 1 , is proportional to with (3 = (3 — a)/2 ~ 0.5. Hence the 
synchrotron emissivity becomes maximum at the cut-off energy and the above results 
suggest the following picture: protons accelerated to ~ i0 20 ~ 21 e y w [\\ De trapped in 
afterglow shock and radiate most of their energy in the TeV range within a time scale of 
a few days. Since the total energy of protons could be comparable to the total kinetic 
energy of a fireball while electrons carry much smaller energy, the energy radiated by proton 
synchrotron around TeV range would be much larger than the energy released as GRB 
photons. 

Now we proceed to estimate of the luminosity of proton synchrotron radiation. Suppose 
that accelerated protons in shocked matter have a power-law spectrum (a = 2) in the range 
7o < % < lu ~ 10 12 and total kinetic energy carried by them in the observer's frame is ^ p E: 

dNp = tjpE 2 

d% m p ln( 7M / 7o 2 ) 7p ' 1 ] 

where 7p = 77p is the proton Lorentz factor at the observer's frame and £ p ~ 1 if the 
accelerated protons are in equipartition. The observed luminosity is given by 

dL 2 d%dN p . 

L( ^"^ =27£ ^rf7; js - p ' (2) 

where j sy n,p = 4aT^ 2 f/mag7p/(3m 2 ) is the synchrotron energy loss rate of a proton in the 
shell frame. After some calculations we find 

Lie,) = 2.5 x 10 45 UTE!l%liW/%!X%y ^ s" 1 , (3) 

where we have assumed ( 7 o,7«) = (10 2 ,10 12 ). If we observe this emission from 
a distance of d = 3000^3 Mpc (z ~ 1), then the observed flux above 1 TeV is 
5.9 x 1CT 10 ^pCb 4 ^^ 4 ^^!^ 2 ^ 3 / 4 ^ 2 photons cm" 2 sec" 1 . This flux is further attenuated 
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by the e creation with intergalactic infrared photon field. The current estimate of the 
optical depth for this intergalactic absorption is still highly uncertain, but a typical value 
for TeV gamma-rays is r ~ 10 (e r = 2.2 x 10 4 ) for z ~ 1 (Salamon & Stecker 1998). The 
amount of infrared background is related to the amount of stars in the universe and this 
is uncertain by a factor of about 2. A factor of 2 reduction of the estimate of r results in 
the attenuation of e T ~ 150. Therefore the above flux would be attenuated by a factor of 
at least 100, and the attenuated flux is consistent with the upper limits set by the Whipple 
telescope (Connaughton et al. 1997) for some GRBs, which are about 10~ 10 -10~ 9 cm~ 2 sec _1 
depending on the source position in the field of view (~ 3°). However, if a burst location 
is determined as well as some of recent GRBs for which afterglows are detected, and 
observation is made for a time scale of day, a flux of ~ 10~ 12 cm _2 sec _1 is detectable 
by currently working ground-based air Cerenkov telescopes (see, e.g., Kifune 1996 for a 
general review), and hence the TeV photons from GRBs at z ~ 1 are marginally detectable. 
Detectability increases rapidly with decreasing distance because of the decrease of optical 
depth as well as increase of the original (unattenuated) flux, and TeV gamma-rays from 
z ~ 0.5 would be easily detectable. 

The above estimate is based on the relatively small distance scale of GRBs, z max ~ 1, 
but larger distance scales are also suggested by cosmic evolution of star formation rate 
(Totani 1997, 1998; Sahu et al. 1997; Wijers et al. 1998) or by the recently detected host 
galaxy for GRB971214 (Kulkarni et al. 1998). The original TeV flux expected on the Earth 
before absorbed in intergalactic fields is almost insensitive to the unknown distance scale of 
GRBs, because we have just scaled the total energy (E) from the observed energy emitted 
as GRB photons. The increase of the intergalactic optical depth with z beyond z ~ 1 is 
also rather slow compared to z < 1 (Salamon & Stecker 1998), and hence the detectability 
of TeV gamma-rays is not so sensitive to the GRB distance scale, if z max > 1. More precise 
estimate of detectability requires better determination of infrared background, and in other 
words, discovery of the TeV afterglow would give important infomation for the intergalactic 
infrared photon field. 



4. Detectability of GeV Photons 

GRB940217 has the third-largest energy fluence in the 4B BATSE catalog (Paciesas 
et al. 1997), and this GRB is famous for the detection of high energy photons by the 
EGRET detector with very long duration (Hurley et al. 1994). The EGRET detected high 
energy photons ranging from 36 MeV to 18 GeV during ~ 5000 seconds. We show that this 
EGRET photons are well explained by the proton synchrotron of the model. The observer's 
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time when the external shock phase begins, t d , is given by J"d/(27g), where the deceleration 
radius is r d = (Y7bE /167mm p 7 2 ) 1//3 . This deceleration time sensitively depends on 7 , and 
it can be as short as t d = 1.3 &2oo-^52 3n ]~ 1//3 7iooo 3 sec when 7 ~ 1000. Therefore the EGRET 
photons can be considered as the external shock origin and our model is applicable, although 
other explanations by internal shocks may also be possible. The observed photon spectrum 
(Fig. 3 of Hurley et al. 1994) seems consistent with the standard spectrum of synchrotron 
radiation, dn/de^ oc £~ 3//2 , and by fitting the data with this photon index we find that the 
differential photon flux dn/de^ is ~ 2 x 10" 11 and 2 x 10 _12 e 7 Qg V photons cm -2 sec -1 keV -1 
for the first 180 sec and delayed photons (180-5400 sec), respectively. This time evolution 
is consistent with the t~ 3//4 profile of Eq. (|3p. If we assume that the fluence of this GRB in 
the BATSE range, 6.6 x 10~ 4 erg cm -2 , is 1/2000 of the total energy E, the distance to this 
GRB is d — 113 &2oo^52 2 Mpc and hence the differential photon flux obtained from Eq. (|3]) 
is 5.6 x 10 _11 ^£ i &200% ^5 £ 7,GeV photons cm -2 sec -1 keV" 1 , where t 5 = t o b s /(5000 
sec). This photon flux is consistent with the observation if the energy conversion into 
accelerated protons and magnetic field is near the equipartition: £ p ~ ^ 4 ~ 0.2nx 1//4 . 
Therefore the delayed GeV photons from GRB940217 are naturally explained by our model. 
On the other hand, there exist some GRBs which are as bright as the GRB940217 but not 
accompanied by such long-duration GeV photons. In such GRBs, the onset of external 
shock phase might be very long after the GRBs and/or the density of interstellar matter is 
quite low. In fact, if the progenitor of a GRB is a massive star, the intensive stellar wind 
prior to the death of the star could have swept up the interstellar medium near the star. In 
this case the TeV or GeV luminosity which is proportional to n 1 ^ 2 could be very small and 
significantly delayed compared to GRBs. 

Because of the e ± pair-creation in intergalactic field, more than 99 % of TeV 
gamma-rays from z ~ 1 must disappear before reaching the Earth. The created e pairs, 
whose energy is about ~ TeV, lose their energy by the inverse- Compton (IC) scattering 
of the cosmic microwave background photons and typical energy of the secondary IC 
photons is 62 ~ 0.6e 2 TeV GeV, which is in the detectable range of the EGRET. The 
expected time delay of these secondary photons is ~ d/(2cy 2 ir ) = 1.8 d^^ 2 day (Cheng & 
Cheng 1996), where 7 pa i r = 10 6 76 is the Lorentz factor of created pair. If the attenuation 
of TeV gamma-rays is significant (e r 3> 1), almost all energy originally emitted in TeV 
range should be converted into GeV range, which is much larger than the original energy 
emitted in GeV photons by proton synchrotron. This effect becomes significant with 
increasing optical depth for TeV photons and compensate the decrease of flux due to 
the increase of distance, and hence we might be able to detect delayed GeV photons 
for rather distant GRBs. In the limit of e T ^> 1 and neglecting the time delay due to 
propagation in the intergalactic field, we have estimated the differential photon flux as 
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1.7 x 10~ 7 £, P £, 3 B E^ 4 b%o n i ^dfy ^3~ 2£ 2,GeV photons GeV -1 cm -2 s -1 . This estimate may be 
further reduced by the delay of ~ a few days, but not so far from the EGRET sensitivity. 
Delayed GeV emission on a time scale of a few days from GRBs at cosmological distances 
may be detectable by the EGRET, or is likely to be detected by the future GLAST 
experiment. Future ground-based telescopes with reduced threshold energy down to tenth 
of GeV range, e.g., the VERITAS project (Weekes et al. 1998) will also be useful for search 
of the secondary GeV photons. 

5. Discussion 

The typical Lorentz factor of electrons in afterglow shock is 7 e = £e(mp/^e)7 and we 
have considered the case of £ e ~ 1/2000. Then the observed synchrotron photon energy of 
a electron with a Lorentz factor 7 e is 77 2 el?/m e = 2.5 x 10~ 4 (2000C, e ) 2 ^bK^El^t^y 2 e V, 
which is the radio band, and it seems to contradict the observations of X-ray or optical 
afterglows. However acceleration of electrons and/or IC scattering of synchrotron photons 
can raise the photon energy. Furthermore, £ e may also increase with time in afterglow. In 
fact, T ep /(r/7) = 870&2oo /3 -E52 1//3n r 2 ^ 3 7 8 ^ 3 decreases as oc 7 s / 3 oc t~ l , and it is possible that 
energy transfer from protons into electrons becomes efficient gradually with the expansion 
of afterglow shock. Note that there are considerable variations for the behavior of afterglows 
observed in X, optical, and radio bands (e.g., Groot et al. 1998). Efficiency of energy 
transfer into electrons and its time evolution could have large variations among GRBs and 
it may be one of the origins of the complicated behavior of GRB afterglows. 

It should be noted that the proton synchrotron emission extends to X-ray, optical, 
and radio bands with the standard synchrotron spectrum of dL/de 1 oc e( 1_Q )/ 2 . If 
a = 2, luminosity per decade of photon energy at e 7 = 1 keV is VlO 9 times smaller 
than that in the TeV range, and the differential flux observed from a distance of 3000 
Mpc is ~ 1.5 x lO-^TEll^iy^X^^Us 2 erg cm" 2 sec" 1 keV^ 1 . This flux is 
comparable to the observed flux of X-ray afterglows for distant bursts such as GRB970402 
or GRB970508 (Piro et al. 1997a,b). Therefore the proton synchrotron radiation could 
contribute to the X-ray afterglows, although it is rather difficult to detect in the optical 
or lower energy bands due to the hardness of the spectrum. Note that optical afterglows 
are associated only to a small fraction of GRBs for which X-ray afterglows are detected. 
Therefore it can be speculated that proton synchrotron was dominant in such GRBs. The 
complicated behavior of afterglows may be a consequence of complicated mixture of proton 
synchrotron and electron synchrotron or inverse Compton scattering. 

We finally note that energy emitted as 10 20-21 eV protons must be roughly the same 
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with that emitted as GRB photons, if the GRB is the origin of ultra high energy cosmic 
rays (UHECRs) observed on the Earth (Waxman 1995; Vietri 1995). On the other hand, 
in our model, energy distributed to such protons is much greater (at least by a factor of 
~ 100) than GRB photons. However, as we have shown, such protons are likely trapped in 
afterglow shock and lose their energy by synchrotron radiation. If the escape fraction of 
protons just cancels the overproduction of 10 20 eV protons, GRBs could still be the origin 
of UHECRs. If the escape fraction is further smaller, then the UHECRs must be explained 
by other sources. The Larmor radius becomes larger with increasing proton energy, and 
the escape fraction may increase with proton energy. This suggests a possibility that the 
spectrum of UHECRs becomes significantly harder above 10 20 eV, which should be tested 
by future experiments. 
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